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Abstract 

Integrated photonic devices are poised to play a key role in a wide variety 
of applications, ranging from optical interconnects [lj and sensors [2] to quan¬ 
tum computing [3]. However, only a small library of semi-analytically designed 
devices are currently known [4]. In this paper, we demonstrate the use of an in¬ 
verse design method that explores the full design space of fabricable devices and 
allows us to design devices with previously unattainable functionality, higher 
performance and robustness, and smaller footprints compared to conventional 
devices [5j. We designed a silicon wavelength demultiplexer that splits 1300 nm 
and 1550 nm light from an input waveguide into two output waveguides, and 
fabricated and characterized several devices. The devices display low insertion 
loss (2 — 4 dB), high contrast (12 — 17 dB), and wide bandwidths (~ 100 nm). 

The device footprint is 2.8 x 2.8 /nn, making this the smallest dielectric wave¬ 
length splitter to date. 

Electronic hardware description languages, such as Verilog and VHDL, are widely 
used in industry to design digital and analog circuits M- The automation of 
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large scale circuit design has enabled the development of modern integrated circuits 
which can contain billions of transistors. Photonic devices, however, are effectively 
designed by hand. The designer selects an overall structure based on analytic theory 
and intuition, and then fine tunes the structure using brute-force parameter sweep 
simulations. Due to the undirected nature of this process, only a few degrees of 
freedom (2 — 6) are available to the designer. The field of integrated photonics would 
be revolutionized if the design of optical devices could be automated to the same 
extent as circuit design. 

We have previously developed an algorithm which can automatically design arbi¬ 
trary linear optical devices [5j. Our method allows the user to design by specification, 
whereby the user simply specifies the desired functionality of the device, and the 
algorithm finds a structure which meets these requirements. In particular, our algo¬ 
rithm searches the full design space of fabricable devices with arbitrary topologies. 
These complex, aperiodic structures can provide previously unattainable function¬ 
ality, or higher performance and smaller footprints than traditional devices, due to 
the greatly expanded design space p si EUDunununtti]. Our algorithm uses 
local-optimization techniques based on convex optimization pQ)j to efficiently search 
this enormous parameter space. 

Here, we demonstrate the capabilities of our inverse design algorithm by designing 
and experimentally demonstrating a compact wavelength demultiplexer on a silicon- 
on-insulator (SOI) platform. One of the key functions of silicon photonics is wave¬ 
length division multiplexing (WDM), which multiplies the data capacity of a single 
optical waveguide or fiber optic cable by the number of wavelength channels used 
nsnmiiE]. Unfortunately, conventional wavelength demultiplexers such as arrayed 
waveguide gratings (19], echellc grating demultiplexers [20], and ring resonator arrays 
EH are fairly large, with dimensions ranging from tens to hundreds of microns ra- 
Our device has a footprint of only 2.8 x 2.8 /irn, which is considerably smaller than 
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any previously demonstrated dielectric wavelength splitter [23]. 

Let us now consider the general formulation of the inverse design problem for 
optical devices. We choose to specify performance of our device by defining the mode 
conversion efficiency between sets of input modes and output modes at several discrete 
frequencies. These modes and frequencies are specified by the user, and kept fixed 
during the optimization process. In the limit of a continuous spectrum of frequencies, 
any linear optical device can be described by the coupling between sets of input and 
output modes, making this a remarkably general formulation [24]. 

Suppose the input modes i — 1... M are at frequencies cjj, and can be represented 
by equivalent current density distributions J,;. Then the electric fields E* generated 
by the input modes should satisfy Maxwell’s equations in the frequency domain, 

V x /^V x E i — c e Ej = — iUiJi, (1) 


where e is the electric permittivity, and fiQ is the magnetic permeability of free space. 

We can then specify N t output modes of interest for each input mode i. We define 
the output mode electric fields £ VJ over output surfaces <%, where j = 1... iV,;- The 
device performance is then specified by constraining the amplitude coupled into each 
output mode to be between a VJ and This leads to the constraint, 


a ij — 



4 • E,;d5* 


— Pij 


( 2 ) 


where we have used overlap integrals to compute the mode coupling efficiency into 
each output mode, and assumed that the input and output modes are appropriately 
normalized. 

The inverse design problem thus reduces to finding the permittivity e and electric 
fields E, which simultaneously satisfy physics, described by equation[lj and the device 
performance constraints, described by equation [2] In general, we also have additional 
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constraints on the permittivity e due to fabrication limitations. 

We use two methods for solving this problem, the objective first method [5] and a 
steepest descent method. In the objective first method, we constrain the electric fields 
Ej to satisfy our performance constraints in equation [2j but allow Maxwell’s equations 
to be violated. We then minimize the violation of physics using the Alternating 
Directions Method of Multipliers (ADMM) optimization algorithm [5], We call this 
method “objective first” since we are forcing the fields to satisfy the performance 
objectives first, and then attempting to satisfy Maxwell’s equations. 

In our steepest descent method, we constrain our electric fields Ej to satisfy 
Maxwell’s equations, and define a performance metric function based on the violation 
of our device performance constraints in equation [2] We then compute the local gra¬ 
dient of the performance metric by solving an adjoint electromagnetic problem, and 
perform steepest-gradient descent optimization mm- 

To design the compact wavelength demultiplexer, we chose a simple planar 3-port 
structure with one input waveguide, two output waveguides, and a square design re¬ 
gion, as illustrated in figure [l]a. For ease of fabrication, the structure was constrained 
to a single fully etched 220 nm thick Si layer on a SiC >2 substrate with air cladding. 
Refractive indices of nsi = 3.49, nsio 2 = 1.45, and n_\\ r = 1 were used. The funda¬ 
mental TE mode of the input waveguide was used as the input mode for the inverse 
design procedure, and the fundamental TE modes of the two output waveguides were 
used as the output modes. At 1300 nm, we specified that > 90% of the input power 
should be coupled out of port 2 and < 1% of the power should be coupled out of port 
3; the converse was specified for 1550 nm. 

The optimization processes proceeded in several stages, as outlined in figure 0>- 
In the first stage, the permittivity e in the design region was allowed to vary con¬ 
tinuously between the permittivity of silicon and air (linear parameterization). The 
objective first method was used to generate an initial guess for the structure, and 
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then the steepest descent method was to fine time the structure. In the second stage, 
the structure was converted to a binary level-set representation [25], and then opti¬ 
mized using steepest descent. Up to this point, the device performance had only been 
specified at the two center wavelengths, 1300 nm and 1550 nm. In the final stage, 
the device was optimized for broadband performance by specifying the device per¬ 
formance at 10 different wavelengths, with 5 frequencies equally spaced about each 
center frequency. Broadband performance was previously shown to be a heuristic 
for structures which are tolerant to fabrication imperfections, and it was hoped that 
this would result in a more robust design [5]. The WDM device was designed in 
approximately 36 hours using a single server with three NVidia GTX Titan graphics 
cards. 

The final designed device is shown in [2^. The simulated electric fields at the 
central operating wavelengths of 1300 nm and 1550 nm are plotted in figure [2}o. At 
both wavelengths, the light takes a relatively confined path through the structure, 
despite the convoluted geometry of the etched silicon layer. 

The devices were fabricated by using electron beam lithography followed by plasma 
etching. Scanning electron microscopy (SEM) images of a final fabricated device is 
shown in figure [3^. The original design was accurately reproduced by the fabrication 
process, with the exception of two small (~ 100 nm) holes next to the input waveguide 
which are missing. 

The measured and simulated S-parameters for the compact WDM device are plot¬ 
ted in figure [4] The plotted wavelength range was limited by the spectral bandwidth 
of the excitation source. Measurements from three identically fabricated devices are 
plotted together in figure it- showing that device performance is highly repeatable. 
Although somewhat degraded in performance with respect to the simulated devices, 
the fabricated WDM devices exhibit relatively low insertion loss (2 —4dB), high 
contrast (12 — 17 dB), and very broadband (~ 100 nm) pass and stop bands. We 
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attribute the discrepancy between simulation and measurement to fabrication imper¬ 
fections. 

In summary, we have experimentally demonstrated a compact, practical wave¬ 
length demultiplexer designed using our inverse design algorithm. This device pro¬ 
vides functionality which never before been demonstrated in such a small structure. 
Our results suggest that the inverse design of optical devices will revolutionize inte¬ 
grated photonics, ushering in a new generation of highly compact optical devices with 
novel functionality and high efficiencies. 

Methods 

Optimization Algorithm and Electromagnetic Simulations 

The detailed inverse design algorithm has been previously described elsewhere Ena, 

1201. A graphical processing unit (GPU) accelerated implementation of the MaxwcllFDFD 
finite-difference frequency-domain solver was used to efficiently solve Maxwell’s equa¬ 
tions throughout the optimization process BZIEBI. 

Fabrication 

The devices were fabricated using Unibond™ SmartCut™ silicon-on-insulator (SOI) 
wafers obtained from SOITEC, with a nominal 220 nm device layer and 3.0 /jrn 
BOX layer. A JEOL JBX-6300FS electron beam lithography system was used to 
pattern a 330 nm ZEP-520A electron beam resist layer spun on the samples. We 
did not compensate for the proximity effect in the electron beam lithography step. 

A transformer-coupled plasma (TCP) etcher was used to transfer the mask to the 
silicon device layer, using a C 2 F 6 breakthrough step and a BCI3/CI2/O2 chemistry 
main etch. The mask was stripped by soaking in Microposit Remover 1165, followed 
by a piranha clean using a 4 : 1 ratio of concentrated sulfuric acid and 30% hydrogen 
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peroxide. Finally, the samples were diced and polished to expose the waveguide facets 
for edge coupling. Detailed schematics of the device are available in the supplementary 
information. 

Measurement 

The devices were measured by edge-coupling the input and output waveguides to 
lensed fibers. A polarization maintaining (PM) lensed fiber was used on the input 
side to ensure that only the fundamental TE waveguide mode was excited. The po¬ 
larization extinction ratio of the light emitted by the PM lensed fiber was measured 
using a polarizing beamsplitter to be 19.0 dB at 1470 nm, and 20.7 dB at 1570 nm. 
A non-polarization maintaining lensed fiber was used to collect light from the out¬ 
puts. The lensed fibers were aligned by optimizing the transmission of a laser at 
1470 nm, ensuring consistent coupling regardless of the transmission characteristics 
of the devices. 

A fiber-coupled broadband light-emitting diode (LED) source and fiber-coupled 
optical spectrum analyzer (OSA) were used to characterize the devices. The transmis¬ 
sion measured through each device was normalized with respect to a straight-through 
waveguide running parallel to each device. This eliminated any coupling and waveg¬ 
uide losses, and yielded a direct measurement of the device efficiencies. 
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Figure 1: Overview of the inverse design process, (a) The device functionality is 
defined for the inverse design algorithm by specifying the surrounding structure, the 
design region, and the coupling between a set of input and output modes. For the 
compact wavelength demultiplexer demonstrated in this work, the structure consists 
of one input waveguide, two output waveguides, and a 2.8 x 2.8 /mi design region. 
1300 nm band light is coupled into the fundamental TE mode of port 2, and 1550 nm 
band light is coupled into the fundamental TE mode of port 3. All three waveguides 
are identical, with a width of 500 nm. (b) Intermediate structures generated by the 
inverse design process. In the first stage, the structure is optimized while allowing the 
permittivity e to continuously vary (linear parameterization). In the next stage, we 
convert to a boundary parameterization and optimize the structure for operation at 
only 1300 nm and 1550 nm. In the final stage, we perform a broadband optimization 
to generate a robust device. 
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Figure 2: The compact wavelength demultiplexer designed by the inverse design al¬ 
gorithm. (a) A three-dimensional rendering of the structure. Silicon is coloured grey, 
and SiOo is coloured blue, (b) Field plots of the device operating at 1300 nm and 
1550 nm. Here, we have plotted the electric energy density U = e|A| 2 . 



Figure 3: Scanning electron microscopy (SEM) images of the fabricated wavelength 
demultiplexer. The device was fabricated by fully etching the 220 nm thick device 
layer of a silicon-on-insulator (SOI) substrate, (a) Top down view, (b) Angled view. 
The vertical sidewalls are clearly visible in this view. 
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Figure 4: S-parameters for the device. Here, we have plotted transmission from input 
port 1 to output ports 2 and 3. (a) S-parameters simulated using finite-difference 
frequency-domain (FDFD) simulations, (b) Measured S-parameters for 3 identical 
devices. The shaded areas indicate the minimum and maximum measured values 
across all measured devices, and solid lines indicate the average values. The insertion 
losses and contrast of the device are somewhat degraded with respect to the simulated 
values due to fabrication imperfections. 
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Figure 5: Detailed schematic of the compact wavelength demultiplexer. 



15 




